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Abstract: We simulate the arrival distribution of ultra-high-energy (UHE) protons by following their 
propagation processes in several strengths of a structured extragalactic magnetic field (EGMF). Compar- 
ing our result to observational one by Akeno Giant Air Shower Array, we constrain the number density 
of UHE cosmic ray sources with the small-scale anisotropy. As a result, the source number density is 
~ 10 -5 Mpc -3 with uncertainty of about an order of magnitude due to the small number of observed 
events. This hardly depends on our structured EGMF strength. We also investigate future prospects 
for this approarch. The near future observations, such as Pierre Auger Observatory, can distinguish 
1(T 6 Mpc -3 accurately from the more source densities. Number of events to discriminate between 
1CT 4 and 10~ 5 Mpc -3 is dependent on the EGMF strength. 
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Introduction 

The origin of ultra-high-energy cosmic rays (UHE- 
CRs) is one of challenging problems in astropar- 
ticle physics. One of significant information on 
UHECR sources is their arrival distribution. Akeno 
Giant Air Shower Array (AGASA) reports small- 
scale anisotropy (SSA) within a few degree scale 
while large-scale isotropy (LSI) based on a har- 
monic analysis [1]. 

The SSA is predicted if UHECRs are of astrophys- 
ical origin with very small number. Therefore, the 
SSA has constrained their source number density 
to 10~ 5 Mpc for no extragalactic magnetic field 
(EGMF) [4, 5] and 10" 6 Mpc for a simply uniform 
turbulent EGMF[6]. 

Recent simulations of cosmological structure for- 
mation predict structured magnetic fields [7, 8] 
which roughly trace the baryon density distribu- 
tion. Tuhs, the EGMF is also structured. In last 
year, we discussed the propagation of UHE pro- 
tons in a structured EGMF that reproduces the ob- 
served local structure [9], and the arrival distri- 
bution. However, we discussed only one EGMF 
strength, which is normalized to 0.4/zG at the cen- 



ter of the Virgo cluster. Observations of mag- 
netic fields in a cluster have large uncertainty in 
the range of 0.1- a few pG [10]. Thus, it is very 
important to investigate the propagation and con- 
straints on UHECR sources in several strengths of 
the EGMF. 

In this study, we discuss the propagation of UHE 
protons in several strengths of the EGMF and a 
Galactic magnetic field (GMF), and compare the 
resulting arrival distributions with observational 
results by AGASA. From the comparison, the 
number density of UHECR sources is constrained. 
Such constraint has large uncertainty due to the the 
small number of observed events at present. So, 
we also discuss the possibility of a decrease in the 
uncertainty with future observations. 

Numerical Methods & Model 

The propagation of UHE protons is calculated by 
an application of the backtracking method, which 



1 . Such AGASA results conflict with High Resolution 
Fly's eye(HiRes) reports, which finds no significant SSA 
[2]. However, this discrepancy is not statistically signif- 
icant due to the small number of observed event [3]. 
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is a method developed in our paper [9]. It is very 
insufficient to calculate their propagation forward 
since cosmic rays do not always reach the Earth 
under finite EGMF even if they are injected toward 
the Earth from their sources. 

Our models of the source distribution and a struc- 
tured EGMF are constructed out of the Infrared As- 
tronomical Satellite Point Source Catalogue Red- 
shift Smvey (IRAS PSCz) catalog of galaxies[ll]. 
This catalog has very large sky coverage (about 
84% of all the sky). Thus, our source distribution 
and EGMF structure reflect large scale structures 
actually observed within 100 Mpc. The EGMF 
strength is normalized at the center of the Virgo 
cluster to 0.0,0.1,0.4 and 1.0 ^G. Outside 100 
Mpc, EGMF is assumed to be an uniform tur- 
bulence with InG and the source distribution is 
isotropic. More details are written in reference [9]. 

In this study, we adopt only a source model that all 
sources have the same power. In conclusion, we 
discuss results from another simple source model, 
which power of each source is proportional to its 
luminosity. 
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Figure 1: xi s as a function of the source number 
density. The error bars originate from 100 times 
source selection. The GMF is considered in the 
lower panel while not in the upper panel. 



Results 

As written above, the arrival distribution has im- 
portant information on UHECR sources. We in- 
vestigate the number density of UHECR sources 
which can best reproduce the AGASA results. Our 
calculated arrival distributions are compared to the 
observed arrival distribution with the two-point 
correlation function 

N V) = l —r Q i a za i Y, 
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which is an indicator of SSA. Number of cosmic 
ray events is set to 49 events in the energy range of 
4 x 10 19 < E < 10 20 cV. For the comparison, we 
define X e max as 
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where N{8) is the two-point correlation function 
and a (8) is lex error due to the finite number of 
events. Small xio provides good agreement with 
the observational result. 



Xios are shown in figure 1. While the source 
number density with 10~ 7 Mpc -3 results in larger 
value, the others are consistent with each other 
within la statistical error. Hence, only the SSA 
cannot constrain the source number density suffi- 
ciently. 

The arrival distribution must also satisfy the LSI. 
We calculate the two-point correlation function 
again, but from merely source distribution to be 
able to predict LSI observed by AGASA. This is 
figure 2. 

In figure 2, the middle panels are the number den- 
sities that best reproduce AGASA results. How- 
ever, source number densities an order of magni- 
tude more than those of the best fit are consistent 
with the observation within la error except B = 
O.O^tG. On the other hand, almost all of source 
distributions with 10~ 7 Mpc~ 3 cannot satisfy the 
LSI. For B = 1.0/xG, there is no source distribu- 
tion in 100 source distributions. This fact can be 
understood in figure 1. Thus, the source number 
density that can best reproduce the AGASA result 
is 10~ 4 - IO" 5 Mpc~ 3 for B = 0.0, O.l^G, and 
IO" 5 - 10~ 6 Mpc~ 3 for B = 0.4,1.0/zG with 
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Figure 2: The two-point correlation functions calculated from only source distributions that can predict the 
LSI. The histograms are the AGASA result within 4 x 10 19 < E < 10 20 eV (49 events). The error bars 
are from the event selection for the finite events and the shaded regions show total la statistical errors. The 
GMF is included. 



uncertainty of about one order of magnitude. The 
source density hardly depends on EGMF strength 
since 95% of space within 100 Mpc has not mag- 
netic field. 

The SSA and the LSI enable us to constrain the 
source number density. However, it has large un- 
certainty which originates probably from the small 
number of observed events. Therefore, one of our 
next interests is how small the uncertainty becomes 
at the Auger era. 

In order to investigate this, it is necessary to com- 
pare our arrival distributions with future observa- 
tional results, which, of course, cannot be known. 
In this study, an isotropic arrival distribution is 
adopted as a template for the future results. If 
UHECR sources are of astrophysical origin and 
have a small number density, the SSA becomes 
stronger. In this viewpoint, we compare our results 
of the simulation to an isotropic distribution, using 
the two-point correlation function. 



Figure. 3 shows distributions of x 2 defined as 

2 = J_ ° v^ aX [N sim (9) - N iso (0)} 2 m 
~ ft 2-^i m\ 2 , to\ 2 ' ( ' 

This value represents the goodness of the fitting. 
The upper panels show current status correspond- 
ing to the AGASA result. The three distributions 
with 10~ 4 , 10~ 5 and 10~ 6 Mpc" 3 are almost over- 
laped. The determination of the number density 
has large uncertainty. 

200 event observation allows us to discriminate 
10~ 5 and 10" 6 Mpc~ 3 with good accuracy since 
it can separate the distributions. This event num- 
ber is probably comparable with current status of 
Auger. Detection of more events can divorce dis- 
tributions with 10~ 4 and 10" 5 Mpc~ 3 . When 500 
events are observed, the number densities of 10~ 4 
and 10~ 5 Mpc" 3 are perfectly separated if EGMF 
does not exist or is very weak. 
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Figure 3: Distributions of x 2 s, calculated from the arrival protons above 4 x 10 19 eV, at several strengths 
of the EGMF. The GMF is considered. The strengths of the EGMF are O.O(left), 0.1, (middle), and 
l.0fj,G(righf). The numbers of events are set to be 49 events within —10° < S < 80° (upper), 200 events 
(middle), and 500 events (lower), within the southern hemisphere to emulate Auger. 



Conclusions 

In this study, we constrain number density of 
UHECR sources with the SSA in several strengths 
of EGMF and investigate future prospects for this 
approach. At the AGASA era, the source num- 
ber density is <~ 10~ 5 Mpc~ 3 with uncertainty of 
about an order of magnitude due to small number 
of observed events. 

That near future observations increasing observed 
event number can improve its uncertainty. 200 
event observation above 4 x 10 19 eV can distin- 
guish 10~ 6 Mpc~ 3 from the more source density. 
This event number is consistent with the number 
observed by Auger until this summer! More event 
detection enables us to estimate the source num- 
ber more accurately and, then, to be easy to com- 
pare it with that of known powerful objects. Num- 
ber of events to discriminate between 10 -4 and 
10~ 5 Mpc~ 3 is dependent on the EGMF strength. 

Finally, we discuss results of another source model 
that power of each source is proportional to its lu- 
minosity, as discussed in [9]. The latter model 
predicts 10~ 4 ~ 10~ 5 Mpc at present. The 
source number density increases since dark sources 
are also counted, but hardly contribute the arrival 
cosmic rays. More observation can discriminate 
10~ 3 Mpc~ 3 from less number densities. This 
model has an additional degree of freedom by lu- 
minosity, compared with the former model. This 
provides large dispersion to distribution of \ 2 - 



Therefore, 10 4 and 10 5 Mpc 3 cannot be dis- 
tinguished even at 500 event observations. 
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